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Abstract 

The special nonlinear dynamical regimes, "bushes of normal modes" , can exist in the 
■ TV-particle Hamiltonian systems with discrete symmetry [Physica D 117 (1998) 43]. The 

dimension of the bush can be essentially less than that of the whole mechanical system. 
' One-dimensional bushes represent the similar nonlinear normal modes introduced by 

CN ■ Rosenberg. A given bush can be excited by imposing the appropriate initial conditions, 

and the energy of the initial excitation turns out to be trapped in this bush. 

In the present paper, we consider all possible vibrational bushes in the simple 
octahedral mechanical system and discuss their stability under assumption that the 
interactions between particles are described by the Lennard-Jones potential. 

s 

: 1 Introduction 

o ■ 

^ ■ A new concept, "bushes of normal modes", was introduced for nonlinear mechanical systems 
with discrete symmetry in fj]. A given bush represents a certain superposition of the modes 
\ associated with different irreducible representations (irreps) of the symmetry group G of the 

& \ mechanical system in equilibrium. The coefficients of this superposition are time-dependent 
functions for which the exact ordinary differential equations can be obtained. In this sense, 
the bush can be considered as a dynamical object whose dimensionality is generally less than 
that of the original mechanical system. The following propositions were justified in previous 
papers Q, |, g|: 

1. A certain subgroup Go of the symmetry group G corresponds to a given bush, and 
this bush can be excited by imposing the appropriate initial conditions with the above 
symmetry group Go C G. 

2. Each mode belonging to the bush possesses its own symmetry group which is greater 
than or equal to the group Gd of the whole bush. 

3. In spite of evolving mode amplitudes, the complete collection of modes in the given 
bush is preserved in time and, in this sense, the bush can be considered as a geometrical 
object. 



4. 



The energy of the initial excitation is trapped in the bush, i.e. it cannot spread to the 
modes which do not belong to the bush, because of the symmetry restrictions. 



5. As an indivisible nonlinear object, the bush exists because of force interactions between 
the modes contained in it. 

6. Taking into account the concrete type of interactions between particles of the consid- 
ered mechanical system can only reduce the dimension of the given bush. 

7. The extension of the bush can be realized as a result of the loss of its stability which 
is accompanied by spontaneous breaking of the bush symmetry (dynamical analog of 
phase transition). 

The special group-theoretical methods for finding bushes of modes are discussed in 0, || . 
The computer implementation of these methods || (see also 0) allowed us to find bushes 
of modes for wide classes of mechanical systems with discrete symmetry. In particular, "ir- 
reducible" bushes of vibrational modes and symmetry determined similar non-linear normal 
modes for all N-particle mechanical systems with the symmetry of any of the 230 space 
groups were found in [|J. The bushes of vibrational modes of small dimensions were found 
and classified into universality classes for all mechanical systems with point groups of crys- 
tallographic symmetry in pj. 

Unlike discussions of existence of bushes of modes, a few papers were devoted to the 
problem of their stability (as an example, we can refer to the investigation of bushes of 
vibrational modes for FPU-a chains in [|J). Nevertheless, the problem of stability is one of 
the most important for bush theory. Indeed, we must keep in mind that bushes of modes 
can be considered as new physical objects only in the case where they are stable in the finite 
domains of pertinent parameters of the appropriate dynamical systems. 

The problem of existence of bushes of modes can be studied with the aid of the group- 
theoretical methods only, and appropriate results do not depend on the concrete type of 
interactions between particles of our mechanical system. In contrast, studying of bush sta- 
bility depends on these interactions essentially. 

In the present paper, after the outline of the bush theory in Sec. 2, we consider bushes of 
vibrational modes for an octahedral mechanical system of mass particles depicted in Fig. [I], 
with or without a particle in the center of the octahedron. All bushes of modes are presented 
in Sec. 3. The differential equations describing their dynamics are discussed in Sec. 4. The 
stability of the bushes of vibrational modes is investigated in Sec. 5 for the case of interactions 
between particles described by the Lennard- Jones potential. 

2 Bushes of vibrational modes 

A detailed consideration of the theory of bushes of modes with some theorems about their 
structure and with the appropriate group-theoretical methods can be found in M. An ex- 
tensive description of the computational algorithm for finding bushes in arbitrary crystal 
structures is presented in [|]. Let us now give an outline of the bush theory with emphasis 
on the ideas important for our further discussion. 

2.1 Geometrical aspects 

We consider classical Hamiltonian systems of iV mass points moving near a single equilibrium 
state which can be characterized by a certain point or space symmetry group G. Let the 

3 x iV-dimensional vector, 

X(t) = ( Xl (t),x 2 (£),..., XAr (t)), (1) 



describe the displacements Xj(t) of all particles of our mechanical system from their equilib- 
rium positions. (Here we denote by the three-dimensional vector Xj the displacement of the 
i-th particle along the X, Y and Z axes). 

The vector X(i) can be written as a superposition of all basis vectors <p^ of the irreducible 
representations Tj of the above mentioned symmetry group G: 

x(t) = E^)# ( 2 ) 

i,i 

The coefficients ^ (t) of this superposition depend on time t, while the 3 x iV-dimensional 
time-independent vectors (p^ determine the specific patterns of displacements of all particles 
of our mechanical system. 

Note that the basis vectors tpf are often called "symmetry-adapted coordinates" . In par- 
ticular, they can be normal coordinates. Hereafter, the term, "mode", means an arbitrary 
superposition of basis vectors corresponding to a given irrep Tj. As a result of this defini- 
tion, every term fi^(t)ip^ in the right-hand side of Eq.(pj) is also a mode of the irrep Tj. 
Sometimes, for brevity, we will refer to /4* (t) as a mode, but a reader must imagine that 
this time-dependent coefficient is multiplied by the appropriate 3 x iV-dimensional vector 
<Pi to give the mode in the exact sense. We can also speak about vibrational modes because 
only such type of symmetry-adapted (normal) modes are considered in the present paper. 

Every dynamical regime of the considered mechanical system can be described by the 
appropriate time- dependent vector X(£) which determines a definite instantaneous configu- 
ration of the system. On the other hand, each instantaneous configuration possesses a certain 
symmetry group Gd (in particular, this group may be trivial: Gd = 1) which is a subgroup 
of the symmetry group G of the system in equilibrium (Gd Q G). Moreover, we can also 
ascribe a certain symmetry group to each basis vector cp\ and to each mode corresponding 
to a given irrep Tj (remember that a mode is a superposition of such vectors!), because 
the definite instantaneous configurations correspond to them. The group Gd contains all 
symmetry elements of group G whose action does not change this configuration. 

Let us introduce, as it is usual in group theory, the operators g associated with elements 
g of group G (g E G) which act on 3 x N- dimensional vectors X(£).0 All elements g E G for 
which 

gX(t) = X(t) (3) 

form a certain subgroup Gd Q G, and a complete set of the above operators g (Wg E Gd) 
represents the group Gd- 

It can be shown that the symmetry group is preserved in time in the sense that its 
elements cannot disappear during time evolution. Actually, this property is a consequence 
of the principle of determinism in classical mechanics.! Thus, the equation gX(t) = X(£) 
(g E Gd), or, formally, 

G D X(t) = X(t) (4) 

is valid for every time t. As a consequence, we can classify the different dynamical regimes 
in our nonlinear dynamical system, described by the vectors X(t) from Eqs.(|l],|D, with the 
aid of symmetry groups corresponding to them. 



lr The symmetry elements g G G act on the vectors of three-dimensional Euclidean space. 
2 The phenomenon of spontaneous breaking of symmetry of a given dynamical regime will be considered 
in the next section. 



Using Eq.(^j), we can obtain the similar invariance conditions for each individual irrep Tj 
of the group G (see details in 0): 

(Tj | G D )^ = (5) 

Here (Tj | Gd) is a restriction of the irrep Tj to the subgroup Go of the group G, and 
, fj,$) is an invariant vector of Fj (n^ is the dimension of this irrep). 
To find all modes contributing to a given dynamical regime with symmetry group Go, 
i.e., for the vector X(£) from Eq.(Q), we must solve linear algebraic equations (||) for each 
irrep Tj of the group G. As a result of this procedure, the invariant vector \ij for some irreps 
Tj can turn out to be equal to zero. Such irreps do not contribute to the considered dy- 
namical regime. On the other hand, some nonzero invariant vectors /it ■ for multidimensional 
irreps may be of a very specific form because of definite relations between their components 
(for example, certain components can be equal to each other, or differ only by sign). The 
contributions to X(£) from such irreps Tj possess special forms. In Tables [|, |], we find, for 
instance, dynamical regimes for which the contributions from the irrep r 10 of the group Oh 
are of the forms (t)[<p { ™ } - cp 2 W) } and (t)[cp { ^ 0) + tp { 2 W) + <^ 10) ]. 

Table 1: Bushes of vibrational modes and their stability domains for the octahedral structure 
with Lennard- Jones potential 



N 


Point 
group 


Generators 


Dim. 


Bush of modes 


Range of bush 
stability 


1 


o h 


13,14,33 


1 


A*i Ski 


Z? = 0.001 


2 




2,14,25 


2 


Ml Ski + Ml Ski 


R = 0.009 


3 


D' 2d 


13,38 


3 


J 1 ) J 1 ) + ,/ 5 U 5) + ..(8)^(8) 
Ml SKI i Ml SKI + M3 <K3 


R = 0.028 


4 


G^ v 


14,37 


3 


Ml Ski + Ml SKI + M3 SK3 


i? = 0.003 


5 


D' 2h 


4,13,25 


3 


„WJ1) + ,,(5)^(5) , (7) (7) 

Mi Ski + Mi Ski + M3 SV3 


R = 0.055 


6 




13,28 


5 




i2i = 0.001 
R 2 = 0.002 


7 




9,13 


3 




R = 0.002 


8 


Gs v 


9,37 


3 




i? = 0.025 


9 


Dm 


13,33 


2 




R = 0.078 



Actually, Eq.(H) can be considered as a source of certain selection rules for spreading 
excitation from the root mode to a number of other (secondary) modes. Indeed, if a certain 
mode with the symmetry group Gd is excited at the initial instant (we call it the "root" 
mode), this group determines the symmetry of the whole bush. The condition that the 
appropriate dynamical regime X(£) must be invariant under the action of the above group 



Table 2: Bushes of vibrational modes and their stability domains for the octahedral structure 
with Lennard- Jones potential 



N 


Point 
group 


Generators 


Dim. 


Bush of modes 


Range of bush 
stability 


10 


Dot, 


2 4 25 


3 


Pi ri i Pi Sri i P2 Sr2 


/?' = n 006 


11 


D> 2 


4,13 


4 


pi Ski + pi Ski + P3 SP3 + P3 ¥>3 


# = 0.027 


12 




4, 2d 


5 


(1) , (5) (5) , (5) tS), 
fi\ Vi + P1 Vl + /4 ¥>2 + 


# = 0.001 


13 


( - / 2v 


4,37 


4 


Pi + PI SPl + ^3 <K3 + ^3 ^3 


R! = 0.003 


14 


Co 


13 


7 




R' = 0.001 


15 


c s 


28 


8 


SM> + /4V 5 > + 4V" + *4"*4 7, + 


i?i = 0.013 
R' 2 = 0.002 


16 


C' s 


37 


7 


m io) k o) + ') + 4 io vr' 


# = 0.007 


17 


c 3 


9 


4 


M 8) (*>?> + vf + vf) + 


R' = 0.005 


18 




25 


6 




R" = 0.001 



G D leads to Eq.(£|) and then to Eq.(||). If the vector fj,j for a given irrep Tj proves to be a 
zero vector, then there are no modes belonging to this irrep which contribute to X(t), i.e., 
the initial excitation cannot spread from the root mode to the secondary modes associated 
with the irrep Fj. 

Note that basis vectors associated with a given irrep Tj in Eq.(|2|) turn out to be equal 
to zero when this irrep is not contained in the decomposition of the full vibrational irrep F 
into its irreducible parts Tj. This is a source of the additional selection rules which reduce 
the number of possible vibrational modes in the considered bush. Trying every irrep Tj in 
Eq.([|) and analyzing the above mentioned decomposition of the vibrational representation 
T, we obtain the whole bush of modes with the symmetry group Go in the explicit form. 

Actually, we have just outlined the group-theoretical scheme for finding bushes of vibra- 
tional modes for arbitrary mechanical systems. The detailed consideration of this scheme 
can be found in [|3], |J. 

In conclusion, let us rewrite the vector X(£) from Eq.(^) as the sum of contributions Aj 



from all different irreps Tj of the group G. For each irrep Tj with the dimension rij, we 
introduce a "supervector" 

*i (6) 

which represents the complete set of basis vectors ip^ (i = 1, . . . , n 3 ) of this irrep. Remem- 
bering that the invariant vector fj,j has the form, 

H = {$,..., i42h (7) 

we can rewrite Eq.(^[) in the following compact form: 

x(t) = EA = E(^WI^ / 

j j 

Here we introduced a formal "scalar" product 



(^(*)l*i)=£/4 ,) (t)¥>. W 

i=l 



of n^-dimensional vectors /Lt^ and «f»j (actually, this product represents a certain 3 x N- 
dimensional vector 



2.2 Dynamical aspects 

Let us return to Eq.(Q). We speak about geometrical aspects of the bush theory when con- 
centrating our attention on basis vectors tfi \ an d we speak about dynamical aspects of 
this theory when we focus on time-dependent coefficients [/!f\t), which will be also called 
"modes" (see comments on this term in Sec. [2.1|) . 

If interactions between the particles of our mechanical system are known, exact dynamical 
equations describing the time evolution of a given bush can be written. We will consider this 
question in Sec.[4| of the present paper, while some general statements about the theory of 
bush dynamics are considered here. 

Two types of interactions between modes in nonlinear Hamiltonian system are discussed 
in ||, namely, force interactions and parametric interactions. We can illustrate the difference 
between these types of modal interactions using a simple example. 

Let us consider two different linear oscillators whose coupling is described by only one 
anharmonic term, U = — 7/if/^2, in the potential energy. Dynamical equations for this system 
can be written as follows: 

jli + ulfi! = 2y (Ufa, (9) 

^2 + ^2/^2 = 7A*1- ( 10 ) 

Here 7 is an arbitrary constant characterizing the strength of the interaction of the oscillators. 
We can suppose that Eqs.(PJlOD describe the dynamics of two modes ii\{t) and ^(i) in a 
certain mechanical system. 

An essential disparity between modes (J.i(t) and ^{t) can be seen from the above equa- 
tions. Indeed, if we excite the mode at the initial instant (/ii(to) 0), the mode ^(t) 
cannot be equal to zero (even if it was zero at t = t \) because a nonzero force = 7/^1 
appears in the right-hand side of Eq.flTUp since fii(t) ^ 0. In other words, the dynamical 
regime (J*i(t) ^ 0,/U 2 (£) = cannot exist because of the contradiction with Eq.([H]). Unlike 



this, the dynamical regime ^(t) ^ 0,/ii(t) =0 can exist because such a condition does not 
contradict equations (|9[) and fllCf). We can say that now there is no force in the right-hand 
side of Eq.([10f) because ^ = as a consequence of the identity H\(t) = 0. 

In the last case, the dynamical regime of the system (§,|n3) represents a harmonic oscilla- 
tion only of the second variable: 

/i 2 (t) = Acos(uj 2 t + 5), (11) 

where A and 5 are two arbitrary constants. 

Thus, there is force action from the mode fii(t) on the mode init), but not vice versa. 
We proved in [|J that such a situation can be realized only in the case where the symmetry 
group of the mode Hx{t) is less than or equal to that of the mode /^(^-^ 

Nevertheless, the mode ^(t) can excite the mode fii(t) under certain circumstances. 
Indeed, substituting the solution (|TT1) of Eq. ([l0|) into Eq.(||), we obtain 



/ii(t) + oj{ - 2A^cos{uj 2 t + 5) n x = 0. (12) 

By means of simple algebraic transformations this equation can be converted to the Mathieu 
equation in its standard form: 

z" + [a-2qcos{2r)]z = 0, (13) 



where z = z{r). But in the (a — q) plane of the Mathieu equation (|T3"D, there exist domains 
of stable and unstable movement. If pertinent parameters 0*2,7) °f our dynamical sys- 
tem (PIJTOD have values such that corresponding parameters a and q of Eq.flTB"!) get into an 
unstable domain, then the nonzero function z{r) and, therefore, the mode /ii(t) appears. 
In other words, the initial dynamical regime /Ui(i) = 0,//2(£) ^ loses its stability, and a 
new dynamical regime /Ui(£) ^ 0,/i2(£) ^ arises spontaneously for definite values of the 
parameters of Eqs.(P,[l0|). Since this phenomenon is similar to the well-known parametric 
resonance, we can speak, in such a case, about parametric action from the mode ^(t) on 
the mode Hi(t). 

The characteristic property of the parametric interaction is that the appropriate force 
(—■§j^ = 27^1^2, in our case) vanishes when the mode (pi, in our case), on which this force 
acts, becomes zero. The following important result was proven in 0: the mode of lower 
symmetry acts on the mode of higher symmetry by force interaction, while the mode of 
higher symmetry can act on the mode of lower symmetry only parametrically. Consequently, 
if the parametric excitation of a certain mode does take place, this phenomenon must by 
necessity be accompanied by spontaneous breaking of symmetry of the mechanical system 
vibrational state. 

Thus, the initially excited dynamical regime can lose its stability because of parametric 
interactions with some zero modes and, as a result, can transform spontaneously into another 
dynamical regime, described by a greater number of dynamical variables, with appropriate 
lowering of symmetry. Obviously, we may treat such phenomenon as a dynamical analog of 
a phase transition. 

The above mechanism of loss of stability of bushes of vibrational modes in octahedral 
structures with Lennard- Jones interactions will be discussed in Sec.[5| 

3 Remember, that speaking about a mode we must take into account that our time-dependent coefficient 
is multiplied by the appropriate basis vector of a certain irrep of the group G, and namely this vector 
determines the symmetry group of the considered mode. 



3 Bushes of vibrational modes for the octahedral struc- 
ture 



In this section, we consider bushes of vibrational modes for the octahedral structure without 
a particle in the center of the octahedron depicted in Fig. [I]. This mechanical structure, in 
its equilibrium state, is described by the point symmetry group Oh, whose brief description 
and explicit form of the appropriate irreducible representations can be found in Appendix 
1. The group Oh has ten irreps Tj (j = 1, . . . , 10): four one- dimensional (Ti, T 2 , T 3 , r 4 ), two 
two-dimensional (r 5 ,r 6 ) and four three-dimensional. Note that basis vectors (H), as well 
as invariant vectors (0), depend on the explicit form of the matrix irreps Tj of the group 
G = Oh$ but the contributions Aj (^j) to the vector X(t) from individual Tj do not depend 
on the concrete form of these irreps. 

3.1 Basis vectors of the irreps of the group Oh 

Not all irreps of the group G are contained in the decomposition of the full mechanical 
representation r mech of our system into its irreducible parts. It can be shown that the irreps 
T 2 , T 3 , T 4 and T 6 are not contained in r mech , while T 10 enters into this reducible representation 
twice. One of two copies of T 10 in the decomposition of r mech describes the translational 
movement of the mechanical structure as a rigid body. Considering only vibrational regimes, 
we will exclude this copy of Tio from consideration. Because of the same reason we also 
exclude the irrep Tg which describes rotation of the considered octahedral structured Thus, 
there are only 3 x 6 — 6 = 12 pure vibrational degrees of freedom and, as a consequence, we 
can restrict ourselves by consideration of the 12-dimensional vibrational representation r v ;b r 
instead of the 18-dimensional mechanical representation r mec h. 

Since decomposition of r vibr contains only one copy of each irrep Tj for j = 1,5, 7,8, 10, 
all symmetry adapted coordinates, generated by the basis vectors of these irreps, turn out 
to be the normal coordinates, and we may not distinguish these two notations in the further 
discussion. 

We give the basis vectors of the irreps Tj, which form a basis of the vibrational repre- 
sentation r v ib r , in Table [3[ These basis vectors were found by the well-known method of 
projection operators with the aid of the appropriate computer program. Each of these 18- 
dimensional vectors determines the specific set of displacements (ri| r 2 | r 3 | r 4 | r 5 | r 6 ) of the 
six particles of our octahedral structure (the numbers of particles are given in Fig. [I]). Here 
Tj = (xj,y.j,Zj) is the three-dimensional displacement vector of the j-th particle from its 
equilibrium position determined by its projections along three Cartesian axes X, Y, Z. The 
complete set of basis vectors <p^ from Table |] can be used as a basis for the decomposition 
of an arbitrary vibrational regime X(t) (see Eq.(0)). 

3.2 Bushes of vibrational modes 

As was already discussed in previous sections, a certain symmetry group Go C G corresponds 
to a given bush. Therefore, to find all bushes, we must know all subgroups G$ of the parent 
group G = Oh of our mechanical system in its equilibrium state. Moreover, we must take into 

4 Namely because of this reason we give the explicit form of the irreps of the group Oh in Appendix 1 
(remember that a given irrep is determined up to arbitrary unitary transformation). 

5 Thcre exist also rotational-vibrational bushes for which modes of Tg turn to be root modes, but in the 
present paper, we will study only pure vibrational bushes. 



Table 3: Basis vectors of the irreps of the group Oh for the octahedral structure depicted in 
Fig.l 





Basis vectors 


Symmetry groups 






-1,0,0 


0,1,0 


1,0,0 


0,-1,0 


0,0,-1) 


o h 






1 -1,0, c 


1 0,1,0 


1 1,0,0 


1 0,-1, c 


1 0,0,2) 




(5) 


^ ( 0,0,0 | a/3,0,0 0, a/3, | 


V3,0, 


| 0, -a/3,0 1 0,0,0 ) 


D 2h 


(7) 

Vi 


| ( 0,1,0 | 


0,0,0 | 


0,0,1 | 0,0,0 | ( 


),0,-l | 


0,-1,0) 


D' 2h 


(7) 
<p2 


| ( 1,0,0 j 


0,0,-1 | 


0,0,0 | 0,0,1 | 


0,0,0 | 


-1,0,0) 


D' 2h 


(7) 
¥>3 


| ( 0,0,0 | 


0,-1,0 | 


1,0,0 | 0,1,0 | 


-1,0,0 


0,0,0) 


D' 2h 


(8) 

Vi 


I ( 1,0,0 1 


0,0,0 | 


-1,0,0 | 0,0,0 | 


-1,0,0 


1,0,0) 


D' 2d 


(8) 
¥>2 


1(0,-1,0 


0,1,0 | 


0,0,0 | 0,1,0 | 


0,0,0 | 


0,-1,0) 


D' 2d 


(8) 
¥>3 


1(0,0,0 | 


0,0,-1 | 


0,0,1 | 0,0,-1 


0,0,1 


0,0,0) 


D' 2d 


(10) 

<Pi 




-2,0,0 


1 1,0,0 


1-2,0,0 


1 1,0,0 


1 1,0,0) 


Ci v 


(10) 
¥>2 




0,1,0 


0,-2,0 


1 0,1,0 


1 0,-2,0 


1 0,1,0) 


Civ 


(10) 
¥>3 




1 0,0,1 


1 0,0,1 


0,0,1 


0,0,1 


0,0,-2) 


Civ 



account not only different subgroups, but also their different embeddings (orientations) in 
the parent group Oh. For example, there are three ways of embedding of the group Go = C 2v 
(we denote them by primes on the appropriate Schoenflies symbols: C' 2v , C 2v ) - they differ 
from each other by their generators. We give the generators of all subgroups of the group Oh, 
corresponding to different nonlinear vibrational regimes, in the second column of Tables [i].[2|. 
For the above case we see the following generators for three variants of the subgroup C 2v : 

C 2v : /14, h 2 Q] C' 2v : hi3,h 2S ; C 2v '■ hi,h 3 j. (14) 

Using the list of symmetry elements from ||, one can find that in C 2v the two- fold axis hi 
is the coordinate axis Cf and the mirror plane h 2 e is the coordinate plane a x ; in C' 2v , h±3 is 
the diagonal two- fold axis C 2 y , while the mirror plane h 2 $ is the coordinate plane o~ z ; in C 2v , 
hi is the coordinate axis Cf and the mirror plane /137 coincides with the diagonal plane a Sy - 
All subgroups Gg } of the group G = Oh with their different embeddings and the sets of 



appropriate invariant vectors for each irrep of this group can be found in |10[ B 

Let us remember, that for obtaining a given bush of vibrational modes with symmetry 
group Go ^ G, we must solve the linear algebraic equations fl5|) for each of five irreps 
Tj (j = 1, 5, 7, 8, 10) of the group G = Oh- In such a way, we obtain the invariant vector 
individually for every irrep Tj. It is very essential that for many irreps these invariant vectors 
turn out to be zero vectors. For example, for Gd = -D4/1 there are only two irreps to which 
nonzero invariant vectors correspond: the one-dimensional irrep Fx and the two-dimensional 
irrep T 5 . Moreover, the invariant vector for T 5 has a special form: /z 5 = (a, 0), where a is an 
arbitrary constant appearing in the solution of Eq.(|5]). Therefore, we obtain the following 
expression! for the bush B2[Dih], if replacing the arbitrary constants by time- dependent 

6 The invariant vectors of irreps of group Oh were obtained in the paper [jl0| in connection with discussion 
of some problems of the theory of color symmetry. 
7 A reader can see this result in Table [l]. 



functions (t): 

X(*)= A *f ) (t)rf ) +/xf ) (t)^ 8) . (15) 

(Note that the coefficient before the second basis vector ip^ of the irrep is equal to zero!). 
Analogously, for bush B3[D' 2d \, we find that only three irreps Ti, T 5 and T 8 contribute to the 
corresponding dynamical regime: 

X(i) = t£\t)<p<P + /4 5) (t)^ 5) + ^(t)^. (16) 

Again, we see that a very specific invariant vector fi s = (0, 0, a) corresponds to the three- 
dimensional irrep r 8 . 

Comparing Eqs. (|I~5]) and (0), one finds the identical symbols ^i\t) and ^\t) in two dif- 
ferent bushes B2[D ih ] and B3[D 2d ]. There is no connection between these functions, because 
we exploit here, as well as in Tables |1|J^, the following convention: one and the same nota- 
tion [if* (t) corresponds to different functions when it is used in the description of different 
bushes. 

Let us consider Tables [I] and where all bushes of vibrational modes for the octahedral 
mechanical structure are listed. The root modes of bushes are underlined (remember, that 
root mode possesses the minimal symmetry among all other modes contained in the given 
bush)!. 

Similar to phase transition theory, we can speak about dynamical domains for a given 
bush. Indeed, as it can be seen from Tables |l|fj, there is only one bush with the symmetry 
group Gd = D ih , namely, the bush B2[D Ah \. The 4-fold axis of this group, associated with 
the generator hu, is directed along the Z coordinate axis, i.e., Gd = D± h (see geometrical 
sense of the generators from the third column of Table [I] in the handbook |J). Obviously, 
there must exist also bushes with the same symmetry group D^h, but with 4- fold axis oriented 
along X and Y coordinate axis, i.e. B\D%^ and B[Df h ], because of the equal status of the 
X, Y, Z axes in the parent group G = Oh- The collection of these three equivalenM bushes 
is similar to the set of different domains in the phase transition theoryS. In Tables [E],p| we 
include only one bush from the complete collection of the equivalent bushes. 

As it can be seen from Tables |l],[| the root modes corresponding to some bushes represent 
superpositions of several basis vectors of a given multidimensional irrep. For example, we 
find in the fifth column of this table that the root modes for bushes B7[D^\, B8[Cs v ] and 
B9[D 3 d] are simply sums of all three basis vectors of the irreps T$, T w and T 7 , respectively. 
In terms of invariant vectors, this means that in solving Eq.(^) we obtain for the above 
irreps the solutions of the same form: /z 8 = fi 10 = /j, 7 = (a, a, a), where a is an arbitrary 
constant (this constant is different for each of the irreps r 8 ,r 10 and IV). In writing the 
bushes B7[D 3 ], B8[Cs v ] and B9[D 3d ] we must replace these constants by the appropriate 
functions of time /4 (t) , (t) and fj^ (t), respectively. 

Let us now consider the more complicated cases where the choice of root mode for a 
given bush is not uniquely defined. For example, there are two different variants for choosing 
the root mode for bush BQ[C 2v \- either the sum (^i 8 "* + (f^) °f two basis vectors of three- 
dimensional irrep T 8 , or the difference (v ? i 1 °' ) — v4 1 °' ) ) °f two basis vectors of another three- 



8 A more exact discussion of this topic see below. 

9 Two bushes are equivalent to each other, if there exist a symmetry element g of the parent group G 
which transforms one bush to another. 

10 The total number of such domains is equal to the index of the subgroup Gd in the group G, i.e., to the 
ratio of the orders of these two groups: ||G||/|jG£>||. 



dimensional irrep Tio can represent the root mode.0 Indeed, the excitation of each of these 
modes separately or simultaneously leads to the excitation of the same bush SGfC^]. 

Such a possibility of different choices of the root mode is brought about by the fact that 
both combinations of the basis vectors (cp^ + <p 2 ^) and ~ t P < 2°^) possess the same 

symmetry group Gjj = C' 2v , while the basis vectors of other irreps contributing to bush 
B6[C' 2v \ have higher symmetry, namely, <p^ - Oh, - -D4/1, <p*P - D' 2h (these three point 
groups are supergroups with respect to the symmetry group C' 2v of the root modes). From 
this example we see that all modes of a given bush do have symmetry higher than or equal 
to that of the root mode (see Sec. PTT| ). 

Above, we consider the cases where bushes may be excited by the initial excitation of a 
single mode (even if there are different variants of the choice of this single root mode). But 
there exists bushes of another type, namely, the bushes whose excitation is possible only if 
two different modes, i.e. modes belonging to the different irreps, are excited simultaneously. 
For the considered mechanical system there are four bushes of such a type in Table |2|: 

bii[d' 2 ] r 7 ,r 8; Bi2[c 2v ] r 8 ,r lo; Bi3[C 2 ' v } r 7 ,r lo; Bi7[c 3 ] r 8 ,r 10 . (17) 

Here we also point out the pairs of the irreps whose modes must be excited for the excitation 
of the corresponding bushes. 

Let us consider the bush B13 [C 2v ] in more detail. We can write down the symmetry groups 
of all four modes contributing to it (these groups can be found in Table [|): 

^ O h , ¥>! 5) - D 4h , cpP D' 2h , <pi 10) C Av . (18) 
Note that the symmetry group C 2v of the whole bush B13[C 2v ] is not found among the 



symmetry groups fll8D of its individual modes. This fact can be understood as follows. Since 
our bush is the superposition of four modes tp^ ,<p[ , ip 3 and it is clear that its 

symmetry group must be the intersection of the four symmetry groups corresponding to 
these modes: 

O h (h 2 , h 3 , h 5 , his, h 2 5),D 4h (h 2 , h u , h 25 ),D 2h (h 4 , h 13 , h 25 ) , C^ v {hu, h 37 ). (19) 

Here we write the generators of each group in the parentheses next to its Schoenflies symbol 
(these generators can be found in the third column of Table [l] near the appropriate symbol 
of the symmetry group)0. 

The two first groups from ( |19"D — Oh and — are subgroups with respect to of the both 
last groups, D 2h and C± v : all elements of D' 2h and C± v are contained in and therefore also 
in Oh- Indeed, the group D± h {h 2 , hu, ^25) contains four two-fold horizontal axes (rotations 
h 2 , h 3 , hi 3 , hu) and four vertical mirror planes (reflections h 26 , h 27 , h 37 , /i 41 ) intersecting along 
vertical four-fold axis (rotations /14, hu, /115) coinciding with the Z coordinate axis, and it is 
obvious that the generators of both groups D' 2h and C& v are identical with some elements of 
the group .D4/JH 

Let us now compare the elements of the groups D' 2h {h A , h 13 , h 25 ) and C 4l) (/ii 4 , h 37 ). The 
second generator of C iv is contained in D' 2h since h 37 = /i 13 • h 25 , and the first generator of 
D' 2h is contained in C± v since h\ A = h^ (see footnote 0). Then we can write the extended sets 



n The invariant vectors (a,a,0) and (a,-a,0) correspond to these two cases. 
12 Only for group Oh, we use the extended set of generators. 



13 Using the multiplication table from || for the element of the group Oh, we can obtain the generators of 
D' 2h and C^ v from the generators of as follows: = hf 4 , /i 13 = ft-2^i4 ; ^37 — ^13^25- 



of the generators of the groups C± v and D' 2d as follows: C4 V (hu, h±, h 37 ), D' 2d (hi3, ^25, h 4 , h 37 ). 
The intersection of these two groups is a group with their common elements, i.e., the group 
with the generators /14 and h 37 . In our notation, this is the group C^X^, ^37) associated 
with the whole bush 513 [CgJ in Table 0. 

Other three bushes from the list (|T7|) can be considered similarly to the above discussed 
bush 513 [C£J. 

Unlike the above cases, the modes of only one irrep are root modes for bushes 

B10[D 3h ),BU[C? 2 ),B16[C2,B18[C i ], (20) 

but these modes are determined by several independent coefficients /Jf (t). For example, 
the root mode fJ.i\t)(fi^ + ^(f)<P2 corresponds to the bush 510 [D 2 h\ (This mode is 
generated by the invariant vector of the form (a,b), associated with the two-dimensional 
irrep r 5 ). Analogously, the root mode jji'°\t)(ipi' ^ + ipi^) + /4i°\t)(p^ corresponds to 
the bush 516[C^] (this mode is induced by the invariant vector of the form (a,a,b) of the 
three-dimensional irrep r 10 ). 



4 Exact dynamical equations for bushes of vibrational 
modes 

We consider a mechanical system of six mass points (particles) whose interactions are de- 
scribed by a pair isotropic potential u(r) where r is the distance between two particles. We 
suppose that in the equilibrium state these particles form a regular octahedron with edge a 
which is depicted in Fig. [TJ. Let us introduce a Cartesian coordinate system. Four particles 
of the above octahedron lie in the XY plane and form a square with edge a^. Two other 
particles lie on the Z axis and we will speak about the "top particle" and the "bottom 
particle" with respect to the direction of the Z axis. Obviously, the distance between each 
of these two particles and any of the four particles in the XY plane is equal to a®. 

In the equilibrium state, our mechanical system possesses the point symmetry group Oh- 
All possible bushes of vibrational modes were considered for this system in the previous 
section, and we remember that a certain subgroup Go of the group G = Oh corresponds to 
each of these bushes. 

In this section, we consider u(r) as an arbitrary pair isotropic potential, but for studying 
the bush stability in Sec.|5] we suppose that u(r) is the well-known Lennard- Jones potential: 

<r) = 4~2-- 6 . (21) 

Here A and 5 are certain constants characterizing the value of the repulsive and attractive 
potential parts, respectively. 

The potential energy of our system in its vibrational state can be written in the form 

V(X) = J2 (22) 

i,j 



where Tij IS Si distance between the i-th and j'-th particles. The iV-dimensional vector 
X = (xi(t), a?2(t), . . • , xjf{t)) in Eq. already introduced in Sec.0, determines the dis- 



placements of all particles at an arbitrary chosen instant t. Note that in our present case 
N = 18.0 

Since all particles possess identical masses, which we suppose to be equal to unity, the 
dynamical equations of Newton type can be written as follows: 

Xi = ~, i = l,...,N (N=18). (23) 

It is easy to find such a scaling transformation that both constants A and B in the potential 
energy V(X) (see Eqs. (pl|, [22^ ) will be equal to unity (A = 1, B — 1): 



t = rnlA^l 



(24) 



Because of this transformation, there are no arbitrary constants in Eq.(^) and, therefore, 
the results obtained in the present paper are universal with respect to all Lennard- Jones 
potentials. Supposing that the scaling transformation (|24]) is already done, we will not change 
the notation of the time variable t and space variables Xi in the further discussion. 



4.1 Lagrange equations for bush modes 

For obtaining dynamical equations of bushes of modes, we must transfer from the original 
(old) variables Xi(t) (i — 1, 2, . . . , N) to new variables fij(t) (j = 1,2, . . . , N), which describe 
the amplitudes of modes in the decomposition (0) 

N 

X(f) = X>i(*)¥>i. (25) 

3=1 

In contrast to Eq. (0) we use here a single index j for numbering the basis vectors of all irreps 
of the group G. 

When a certain bush with a symmetry group Gd C G is considered, some amplitudes 
fJLj(t) turn out to be equal to zero (fJ>j(t) = 0) and we can suppose without loss of generality 
that index j from Eq.(^5|) runs over only the first n values (j = 1,2, ... ,n). Here n < N 
for all bushes of modes except for the bush with trivial symmetry group Gd = 1. Thus, a 
given bush represents a dynamical system whose dimension is less than that of the original 
mechanical system (n < N) and its time evolution can be described with the aid of n 
generalized variables fij(t) (j = l,...,n). The Lagrange method seems to be the most 
natural and convenient for obtaining the dynamical equations for bushes of modes. Let us 
consider this idea in more detail. 

For a given n-dimensional bush we can rewrite Eq.(^) in the form 

n 

X(t) = J2»j(t)<Py (26) 

3=1 

Since all basis vectors ipj are known, Eq. fl2"S|) provides explicit expression of each of the iV 
old variables Xi(t) (i = 1,...,N) as a function (linear superposition!) of n new variables 



u In this section, we denote by N the total number of degrees of freedom of the considered mechanical 
system (in other sections, N is the number of particles). 



fj>j(t) (j = 1, . . . , n). Therefore, the potential energy V(fJ,i, . . . , /i n ) in the new variables can 
be obtained from the potential energy V(xi, . . . ,xn) in the old variables by substitution 
Xi(t) (i = 1, . . . , N) from Eq.fl26|) to the function V(x\, . . . , xn). Similarly, the kinetic energy 
T(xi, . . . , xjy) — \Yli=i%i can De transformed with the aid of Eq.(|26|) to the following 
quadratic form in terms of velocities fij(t) of the new variables: 

~ \ n n 

T(p, x , . . . , fi n ) = ~Y,12 ^^(Vj I <Pk) ■ ( 27 ) 

Z j=l k=l 

Here ((fj\(p k ) is the scalar product of two real ^-dimensional basis vectors (fj and ip k . We 
can suppose that the basis vectors if; {j — 1, . . . , N) are orthonormal: 

(¥>>*> = (28) 

(Note that basis vectors from Table |3| do satisfy this condition). 
Taking into account Eq.(^), we can rewrite Eq.(|27|) as follows: 



T(/ii, . . . ,/i„) 

Then with the aid of the Lagrange function L - 
a given n-dimensional bush with the symmetry 



Z 3=1 

= T — V, the exact dynamical equations for 
group Gd C G can be obtained: 



d / dL \ dL ^ (1 ) (30) 

Let us emphasize that these dynamical equations of the bush of modes are exact. Actually, 
this fact is brought about by the symmetry-related restriction 

G D X(t) = X(t) (31) 

which reduces the number of degrees of freedom from iV to n. Indeed, namely the condition 
( pn| ) leads to the vanishing of some terms in Eq.(^) and reduces it to Eq. (|26|) . 

4.2 Examples of dynamical equations for bushes of vibrational 
modes 

The exact equations ( |30"|) for bushes of vibrational modes for the octahedral structure with the 
Lennard- Jones potential are rather complicated. Because of this reason, we wrote a special 
MAPLE-program for obtaining these equations and for studying the stability of bushes of 
modes in modal space (see Sec. |573| ). Let us now write down the very clear and compact forms 
of the exact dynamical equations for some bushes in terms of natural geometrical variables. 

We consider the following three bushes from Table [TJ: Bl[Oh], B2[D ih \ and iMfCy, whose 
dimension is equal to 1, 2 and 3, respectively. Note that symmetry groups of the above bushes 
are connected to each other by the following group-subgroup relation: C iv C D ih C Oh- 

The geometrical forms of our mechanical system in the vibrational state, corresponding 
to the these bushes, can be revealed from the superpositions of modes (given in the fifth 
column of Table |1|). But it is easier to understand the geometrical sense of the above bushes 
immediately from the symmetry groups corresponding to them. 



Indeed, the one-dimensional bush 51 [0^] consists of only one ("breathing") mode: the 
appropriate nonlinear dynamical regime X(t) = ^i\t)(p{^ describes evolution of a regular 
octahedron whose edge a = a(t) periodically changes in time. 

The two-dimensional bush B2[D 4h ] describes a dynamical regime with two degrees of 
freedom: X(t) = ^ (t) ip^ + ^(t)^^. The symmetry group G D = D Ah of this bush 
contains the 4-fold axis coinciding with the Z coordinate axis and the mirror plane coinciding 
with the XY plane. This symmetry group restricts essentially the form of the polyhedron 
describing our mechanical system in the vibrational state. Indeed, the presence of the 4-fold 
axis demands that the quadrangle in the XY plane be a square. Because of the same reason, 
the four edges connecting the particles in the XY plane (vertices of the above square) with 
the top particle lying on the Z axis must be of the same length which we denote by b(t). 

Similarly, let the length of the edges connecting the bottom particle on the Z axis with 
any of the 4 particles in the XY plane be denoted by c(t). In our present case of the 
bush B2[D 4h ], b(i) = c(t) for any time t because of the presence of the horizontal mirror 
plane in the group Go = D^. But for the three-dimensional bush 54[C 4?) ], described by 
X(t) = Hi\t)<p^ + fif\t)ipf^ + fi^ \t)cfi^°\ this mirror plane is absent and, therefore, 
b(t) / c(t). 

Let us also introduce two heights, hi(t) and h 2 (t), corresponding to the perpendiculars 
dropped, respectively, from the top and bottom vertices of our polyhedron onto the XY 
plane. Now we can write the dynamical equations of the above bushes in terms of pure 
geometrical variables a(t) , b(t) , c(t) , hi(t) and h 2 (t). 

We choose a(t) and h(t) = h\(i) = h 2 (t) as dynamical variables for describing the two- 
dimensional bush B2[D 4h ] and a(t),hi(t) and h 2 (t) as dynamical variables for describing 
the three-dimensional bush BA[C± V ]. Using these variables we can write down the potential 
energy for our bushes of vibrational modes as follows: 



Bl[O h ] 

B2[D 4h ] 

BA[C iv \ 



V m {a) = \2u(a) + 3w( v / 2a), 

V B2 (a, h) = 4u(a) + 2u(^2a) + 8u (^/h 2 + f ) + u(2h), 
V Bi (a, h u h 2 ) = Au{a) + 2u(V2a) + Au{b) + 4u(c) + u(/ii + h 2 ), 



where b = ^ + [\h x - \h 2 )\ c=yj^+ (\h 2 - \h) 2 . 

Then with the aid of the Lagrange method, we can obtain the following dynamical equa- 
tions for the above bushes of vibrational modes: 



Bl[O h ] : 






a = 


-Au'(a) - 


V2u'(V2a); 


B2[D 4h ] : 






a = 


-2u\a) - 


V2u'(V2a) -2«'(6)f, 


h = 


- 


-u'(2h); 


BA[C 4v ] : 






a = 


-2u'(a) - 


V2u'(V2a) -u'(6)f - 


hi = 


u\bY\ 


^-u'(h 1 + h 2 ), 


h 2 = 


-«'(c)^ 





(32) 



Thus, we obtain the dynamical equations of our bushes of vibrational modes in terms of 
variables with explicit geometrical sense. Each bush describes a certain nonlinear dynamical 
regime corresponding to such a vibrational state of the considered mechanical system, that 



at any fixed time the configuration of this system is represented by a definite polyhedron 
with symmetry group Go of the given bush. 

We can write dynamical equations for the above bushes in terms of vibrational modes 
as well. In spite of the more complicated form, these equations turn out to be more useful 
for the bush theory, since they allow us to decompose the appropriate nonlinear dynamical 
regimes into modes of different importance for the case of small oscillations - root modes 
and secondary modes of different orders 0. As an example, we write below the dynamical 
equations for the bush B4 [Cy in terms of its three modes, /4 (t) = j(t) (root mode), 
^(t) = ^(t) = v(t) (secondary modes): 



-~Uy/2u'(a) + Au'(V2a) + 2v!(h x + h 2 ) + '—^-(2y/2a + 5h x - h 2 ) + 
6 b 



u'(c) 



(2V2a + 5h 2 -hx)), (33) 



v 



c 

■-(2V2u'(a) + 2u'(V2a) - 2u'(h l + h 2 ) + ^-(V2a - 5/ii + h 2 ) + 
6 6 



u c 



(V2a-5h2 + hi)), (34) 



7 = 7 \ u ( b ) — i u ( c ) — : — • ( 35 ) 



Here 



a = V2(r + n + u);b= yj_ (r + fi + vf + (r + fi - 2v - 3 7 ) 2 ; 

c = ^(ro + Jl + vf + (r + \i - 2v + 3 7 ) 2 ; (36) 
hi = r + /I - 2v - 27; h 2 = r + /i - 2v + 27; r = ^. 

Let us stress that all above equations are valid not only for the Lennard- Jones potential 
(|2l|) , but for any pair and isotropic potential u(r), as well. 

Now, we want to demonstrate that Eqs.( |33| , |34| , |35D possess the specific bush structure 
which is discussed in Sec.|2|. Note that there are the definite relations between the symmetry 
groups of the modes /i(t), v(t) and j(t): 



i), anu ^l\iy- 
G D [fi] =O h D G D [v) = D 4h d G D [j) = C iv 



(37) 



(see Tables ^] and [Xp . Therefore, it follows from the general bush theory that the mode 7(£), 
being of minimal symmetry, must provide some forces in the right-hand side (rhs) of the 
Eqs.(|3"3"D and ((W) even if /i(t) = 0, y{t) = and, analogously, the mode v{t) must provide 



a certain force only in the rhs of Eq. (|33|) when fi(t) = and j(t) = 0. These properties can 
be easily seen from Eqs.([3~3H36|). One can also reveal that the mode //(t), whose symmetry 
is higher than those of the modes i/{t) and 7(£), does not produce any forces in the rhs of 
Eqs.(|34"D and (|35|), if v{t) = and ^(t) = 0. Similarly, one can verify that the modes v{t) and 
Hit) do not produce any forces in the rhs of Eq.(|3^D when /i(t) ^ 0, u(t) ^ 0, but ^(t) = 0. 

Thus, we confirmed the typical bush structure of the dynamical equations (|33H35|) in 
terms of the individual modes. 



5 Stability of bushes of vibrational modes 



All possible bushes of vibrational modes for the octahedral structure depicted in Fig. |I| 
were found by means of group-theoretical methods in the previous section, and these results 



are independent of the type of interactions between individual particles of our mechanical 
system. On the other hand, studying the stability of bushes of modes depends essentially 
on the concrete type of interactions in the considered system, and we suppose that they can 
be described by the Lennard- Jones potential (pT|). Then the bush stability can be analyzed 
with the aid of numerical methods. 



5.1 Setting up the problem 

First of all, we must define what we mean speaking about stability of bushes of vibrational 
modes, since the term "stability" is used in very different senses. Let us excite a given bush 
B[Gd] with the aid of the appropriate initial conditions (see below). We can speak about 
the energy Eq of the mechanical system in this initial state. Remember that the given bush 
is a closed dynamical object, i.e. energy Eq can spread only among its own modes. Let us 
decompose the 3 x JV -dimensional vector X(t), describing our bush B[Go] at a certain 
instant t into the complete set of orthonormal basis vectors from Table ||. As a result, we 
obtain a definite collection of nonzero coefficients fJ>^\t) from Eq.@. We can repeat this 
procedure at other times during the evolution of the system, thereby obtaining the time 
evolution of /4 (£)■ During this time evolution, /4 (£) may pass through zero but will never 
be identically equal to zero. This set of modes forms the excited bush B[Gd] of vibrational 
modes. 

Increasing gradually the energy E of the initial excitation, we can detect the loss of 
stability of the bush B[Go}- Indeed, beyond a certain value of Eq, after some transient time 
interval, the complete set of modes with nonzero amplitudes /J^\t) becomes larger than that 
of the bush B[Go]- Obviously, this phenomenon leads to the appearance of a certain new 
bush B[Gd] which includes the old bush B[Go] and whose symmetry Go is lower than the 
symmetry group Go of the bush B[Go] (Go C Go), because all nonzero vibrational modes 
with symmetry higher than or equal to Go are already contained in B[Go\- 

Thus, the loss of stability of a given bush is accompanied by the spontaneous breaking 
of symmetry of the initial excited dynamical regime, described by this bush. We already 
discussed this phenomenon in Sec.|] and concluded that its cause is analogous to that of 
parametric resonance. 

We can also say this in other words. A given bush B[Go] represents a certain dynamical 
regime in the considered mechanical system. Its modes interact with other modes which do 
not belong to B[Go], but these interactions must be of parametric (not force!) type only 
(see Sec. p.2|) . For the appropriate initial conditions we can get into a region of unstable 



movement. As a result, some new modes are excited which were forbidden by the principle 
of determinism of classical mechanics. Then we can speak about the loss of stability of the 
original bush B[Go] and its transformation into a larger bush B[Go] with Go C Go- 
It was found from numerical experiments that the boundaries of stable (unstable) domains 
for bushes of vibrational modes depend not only on the initial energy Eq, but on the way of 
excitation, as well. Because of this reason, we will fix the way of initial excitation.lll 

Let us first discuss the loss of stability of the bushes from Table [I]. The root mode of each 
of these bushes is determined by a single time-dependent coefficient fif\t) whose initial 
value at t — to we will denote by the symbol /i (/i = $\t )). At the initial instant 
t = to, we fix the coordinates of all particles of the mechanical system in such a way that 
their displacements correspond to the appropriate root mode with amplitude jUq, while their 



15 Similar results can be also obtained for other ways of excitation. 



velocities are equal to zero. Namely this choice of initial conditions determines the above 
mentioned way of excitation of a given bush. 

Using these initial conditions we solve numerically the exact dynamical equations of the 
considered mechanical system with 18 degrees of freedom, and analyze the set of nonzero 
modes ^/\t) in the decomposition (|[) of the vector X(i) obtained as a result of this solving. 
Then we gradually increase the value /i and repeat the procedure just described until the 
number of nonzero modes fif\t), at some value //o = R, becomes larger than that of the 
bush B[Gr>]. We will refer to R as the threshold of stability of the given bush. Obviously, 
in such a way we obtain the upper boundary of the first stability region of B[Gd] in the 
one-dimensional space of all possible values /io (0 < /io < 

The threshold values R, with the appropriate numerical accuracy^ for different bushes 
for the octahedral structure without a centered particle can be found in the last column 
of Table [I]. For example, the threshold R = 0.028 corresponds to bush B3[D' 2d ]. It means 
that below this threshold only three nonzero modes can be revealed in the decomposition 
of the appropriate vector X(£) : ^f\t) (root mode) and fjp (t) , fjp (t) (secondary modes), 
while beyond this value of R some new modes /4 J (t) appear in the decomposition (|2|) of the 
vector X(i). We treat this fact as evidence of the loss of stability of bush B3[D' 2d ] and its 
transformation to a new bush of higher dimension. Note that in the present paper, we do 
not determine what this new bush is, because realization of such a goal requires a sufficient 
time due to some computational errors problems. 

We already discussed the bush B6 [C' 2v ] from Table [1] which can be excited either by an 
initial excitation of the root mode + ip 2 of the irrep T s with the amplitude /4 (t), or by 
excitation of the root mode — (p 2 10 ^ of the irrep Ti Q with the amplitude fii' \t). In such 
cases, we find two different threshold values R\ and R 2 corresponding to the two different 
root modes. 

A more complicated computational problem arises when we analyze the domains of sta- 
bility for some bushes of modes from Table Indeed, the root modes of bushes BlQ[D 2 h\, 
B11[D' 2 ], B12[C 2v ), B13[C£J, BU[C' 2 ), B16[C' S ), B17[C 3 ], include two independent functions 
/4 (£)■ For some of the above bushes these functions are associated with basis vectors of one 
and the same irrep; in other cases they are associated with basis vectors of two different irreps. 
Therefore, we must now study the boundary of the stability domain in a two-dimensional 
space of initial values, say, /ii(to) an d /^(to)) of the two appropriate time-dependent variables 
fi^\t). Since studying of the form of such a boundary, by means of numerical experiments, 
requires a lot of computational time, we restrict ourselves to the determination of only the 
radius R! of a circle which with assurance lies in the stability domain. Namely these radii we 
give in the last column of Table p] as certain estimation characteristics of stability domains 
of the appropriate bushes of vibrational modes. 

In Table p], as well as in Table [l], we encounter the case where a given bush can be excited 
by using different choices of the root mode, but now each of these modes is determined by two 
amplitudes /i-^(t). This is the bush £?15[C S ] whose root modes are /4 (^Vi + /4 8 (t)ip 2 
or /4 (^Va + /4 (*)v4 ■ Because of this reason, two different radii R[ and R' 2 are given 
in Table ^ corresponded to the above two variants of the root mode. 

Finally, the root mode of the last bush 518[Cj] in Table ^| depends on three functions 
fii (t), /J, 2 (t), /4 (t) and, therefore, we must study, in this case, the boundary of the stability 



16 Note that we do not study the other possible regions of stability of the given bush B[Gd] in the present 
paper. 

17 We hope that numerical errors do not exceed unity in the last digits of numbers in our tables. 



domain in three-dimensional space. We give for this bush the radius R" of a sphere which 
with certainty belongs to the stability domain. 

5.2 Numerical results on bush stability for octahedral structure 
without centered particle 

Up to this point our consideration was based on group-theoretical methods which are in- 
dependent of the type of interactions between particles in the mechanical system. Now we 
intend to discuss the findings on bush stability using for this purpose the concrete type 
of interactions described by the Lennard- Jones potential already discussed in Sec. (| (see 
Eq. flU)). 

First of all, we must examine the stability of the octahedral structure depicted in Fig.l 
in the equilibrium state. To this end, we found squares of all eigenfrequencies in the har- 
monic approximation with the aid of the MAPLE package and made their classification in 
accordance with the irreps of the group Oh- We obtained the following results: 

w 2 [r x ] = 61.082; c; 2 [r 5 ] = 13.835; 

u 2 [T 7 ] = 31.498; to 2 [T 8 ] = 15.346; (38) 
a/[r 9 ] = 0; lo 2 [T w ] = 46.843; 0. 

Here ^[Tj] is a square of the frequency corresponding to the irrep Tj. 

The irreps 1^,1^3, r4 and Tq are not contained in the decomposition of the full vibra- 
tional representation r mech of the octahedral structure and, as a consequence, there are no 
vibrational frequencies corresponding to them. 

Note, CL> 2 [r 9 ] and the last value of u> 2 [r 10 ] are equal to zero. They correspond, respectively, 
to pure rotation and pure translation of the mechanical system as a whole, and must be 
excluded from the further discussion because we consider vibrational regimes only. § All 
other w 2 [Tj] are positive and this proves that the equilibrium state of our mechanical system 
is stable. 

Let us also write down the value a of the edge of the regular octahedron depicted in Fig.|I| 
in equilibrium: 

a = 1.117. (39) 

Hereafter, we use such dimensionless units that A = 1 and B = 1 in the formula fl2~T|) for the 
Lennard- Jones potential (see Sec. f|). 

Now we consider our results on the stability of dynamical regimes described by bushes of 
vibrational modes. The computational scheme, outlined in the previous section, and fourth 
order Runge-Kutta method were used for appropriate numerical experiments. The threshold 
values R, R', R" characterizing the size of stability domains for all bushes are given in the 
last column of Tables [I], ^|. It can be seen from these tables that the size of stability domains 
differs very considerably for different bushes: 0.001 < R < 0.078,0.001 < R' < 0.027. 

The one-dimensional bush Bl[Oh], consisting of the breathing mode only, possesses the 
minimal "stability reserve": R = 0.001. Obviously, this is a very small value for R, because 
the displacements of particles corresponding to it are approximately equal to 0.0005, while 
the equilibrium edge ao of our octahedron is equal to 1.117 according to Eq. (|39|) . Thus, 
the bush I?1[0/J exhibits very weak stability for the octahedral structure without a centered 



18 Remember that movements of rotation-vibration type were excluded from our consideration in the 
present paper. 



particle. Nevertheless, if a particle in the center of the octahedron depicted in Fig .[I] is present, 
the threshold value R for this bush can be up to 100-1000 times greater, depending on the 
force properties of the centered particle (see Sec. |5.4|) . 

The maximal value R corresponds to the bush B9[D 3( i] (R = 0.078), demonstrating how 
strong can be the difference between the stability properties of dynamical regimes with 
different character of displacements of individual particles (note that using Tables [TJ, 0, |3|, 
we can imagine the character of movement of the mechanical system for any given bush) . 

Our main result can be formulated as follows: all bushes of vibrational modes possess sta- 
bility domains of a finite size for the octahedral structure with the Lennard- Jones potential. 



5.3 Interactions between different modes 

In many physical papers, pair interactions between modes are discussed (compare the pho- 
non-phonon interactions in the presence of anharmonic terms in the appropriate Hamilto- 
nian). In a rigorous sense, such a setting of the problem is incorrect and can be used only as a 
first step of investigation, for example, to obtain some estimations. Indeed, we can formally 
write the Lagrange equations describing dynamics of only two arbitrarily chosen modes, 
but this new two-dimensional dynamical system may not correspond to any real dynamical 
regime in the original mechanical system. 

The importance of the bush theory is brought about, above all, by the following conse- 
quence: choosing a given bush, we can guarantee that the full collection of its modes is closed, 
i.e. that the appropriate dynamical regime including only these modes can really exist in the 
considered mechanical system. As a consequence, the Lagrange equations, obtained for the 
complete set of the modes of a given bush, are exact. 

Nevertheless, studying the interactions between pairs of different normal modes may be 
useful in certain cases as an approximate method for revealing some essential properties of 
nonlinear dynamics of the mechanical system. Let us illustrate this idea using the following 
examples. 

1. We consider a possibility of the parametric excitation of "sleeping" modest!, belonging 
to the irreps r 5 ,r 7 ,r 8 and r 9 , brought about by their pair interactions with the breath- 
ing mode /ii(t). In other words, we must find the threshold values /4^(0) = /j, c beyond 
which each of these previously sleeping modes can appear because of its interaction with 
the "active" mode Hi (t) ^ 0. The following results were obtained by means of numerical 
experiments: 

T 5 [^\t),^\t)-fi c = 0.020], 
T 7 [^\t),$\t)JP(t) -ijl c = 0.173], (40) 

r 8 [^ 8) (t),^ 8) (t),^ 8) (t)-^ c = 0.001], 

r 10 ^ 10) W,^ 10) W,/if \t) ~Hc = 0.333]. 

Note that modes listed in square brackets next to the symbol of the appropriate irrep possess 
the identical threshold value /i c . 

According to Eq. (|40|) , the minimal threshold value /i c corresponds to the modes of the 
irrep r§ (/i c = 0.001), and this value coincides with the threshold of stability R = 0.001 
of the one- dimensional bush i?l[0/i] pointed out in Table [I]. It is essential, that when we 

19 We mean by this term the modes which do not belong to the bush generated by a given root mode 
(fj±\t), in our case) and which were not excited at the initial instant t = to. 



find this value R for .Bl[0/j] by means of the numerical procedure described in Sec. |5.1| , we 
cannot reveal the cause of such a small value of R. Now we can assert that the bush Bl [O/J 
loses its stability because of the interactions with the modes namely of the irrep T§ and, 
therefore, this bush must transform into a certain bush, containing these modes, such as 
B3 [D' 2 d ] , B6 [C' 2v ] , etc., but not into such bushes as B2[D± h \, S4[C 2 «], etc. 

2. Let us discuss the threshold value R = 0.009 of the stability of the bush B2[D ih \ 
(see Table |l|) in terms of pair interactions between different modes. This bush contains two 
modes: the root mode ^^{t) and the secondary mode (t). We can find the threshold of 
the initial value /4^(0) = A* c beyond which the active mode fii(t) excites, by means of its 
parametric action, other modes which were previously sleeping. (There is one and only one 
mode — fJ^\t) — which is brought into the vibrational process by force originating from 
the active mode fi[ (t).) We obtained the following results: 

r 5 [/4 5) (t)-/^ = 0.273], 

r 7 [/4 7) (t),/4 7) (f) -n c = 0.560; fi { P(t) -fi c = 0.403], 
r 8 [/4 8) (£) -fi c = 0.344; $\t) - fi c = 0.344; fif\t) -fi c = 0.455], 
r 10 [/4 10) W - lie = 0.616; fi ( 2 10 \t) -fM c = 0.616; ^°\t) - fi c = 0.577]. 

Note that unlike the case where the breathing mode was active (see Eq. (f40|) ), threshold 
values |U C for the excitation of sleeping modes by the active mode ^ (t) can be different for 
the different modes of one and the same multidimensional irrep. 

A remarkable fact can be revealed in analyzing the above threshold values given in Eq. (f4T|) . 
Indeed, all /jl c , corresponding to the initial value of the root mode fi^\t), are essentially 
greater than the threshold R = 0.009 for the loss of stability of the considered bush B2[D ih \. 
On the other hand, as we already know, the threshold for the parametric excitation of the 
modes associated with the irrep Tg by action from the active mode /4 (t) is a very small 
value: \i c = /4^(0) = 0.001. This fact suggests that the sufficiently weak stability of the bush 
52(^4^] brought about by the parametric excitation of some sleeping modes originated not 
from the root, but from the secondary mode of our bush! Indeed, if we consider dynamics 
of the bush B2[D^h] when /x^(0) are slightly higher than 0.009, the value of its secondary 
mode ^i 1 ^) does reach the threshold value 0.001 which can lead to excitation of the modes 
of the irreps T 8 . 

Thus, we encounter now the case where a given bush loses its stability because of the 
phenomenon similar to the parametric resonance induced by its secondary mode, and this 
property demonstrates, in particular, that a bush is an indivisible dynamical object. 

5.4 Stability of bushes of vibrational modes for octahedral struc- 
ture with the centered particle 

First of all, let us note that all octahedral molecules, known to us at the present time, 
possess an atom in the center of the octahedron (see Fig.|Tj) . This fact suggests that the 
stability of such structures can be greater than that of the mechanical system considered 
in the previous sections. Taking into account this hint, we examined bush stability for the 
mechanical structure depicted in Fig.|l], supposing that the particle in the center of the 
octahedron is described by the Lennard- Jones potential different from that of six peripheral 
particles. 



(41) 



Let us assume that the Lennard- Jones potential of the form 

U{r) = ± ~ £ (42) 

with B > 1 corresponds to the centered particle, while all other particles are described by 
the potential (|42|) with B = 1. Such an assumption provides us a possibility to make the 
attractive part of the potential of the centered particle greater than that of peripheral atoms 
in spite of the same repulsive part. 

We will consider only several bushes of modes for the centered structure to show that the 
size of the stability domain can be increased essentially by means of the appropriate choice 
of the value B. Note for the beginning of the discussion, that all bushes from Tables Q^j are 
relevant bushes for centered structure, as well as for that without the particle in the center 
of the octahedron, if this particle is fixed.0 

Let us discuss some results of the appropriate numerical experiments. The boundary of the 
stability domain for the bush BlfOJ can be increased from the threshold value R = 0.001 
(see Table p up to the value R = 1.01 for the case of the centered structure with B = 5.5. (It 
is interesting that the function R(B) for this bush is not monotonic). The threshold values 
R for other considered bushes also become larger for B = 5.5 in comparison with those for 
B — 1. For example, we obtained R = 0.011 for bush B2\D^ instead of R = 0.009 and 
R = 0.118 for bush i?4[C 4 ^] instead of R = 0.003. Note that the dependence of the boundary 
R of stability domains on the value B is essentially different for different bushes. 

We can conclude that the stability of bushes of vibrational modes for the centered struc- 
ture can be increased, and for some bushes substantially increased, in comparison to that of 
the mechanical system without the centered particle on account of the appropriate choice of 
the constant B in the Lennard- Jones potential (f£^). 



6 Summary 

The specific dynamical objects — bushes of normal modes — whose theory was developed 
in [l], |2|, ||, are discussed in the present paper for octahedral mechanical systems with point 
masses. Being a certain collection of modes, a given bush possesses some symmetry group 
which is a subgroup of the symmetry group of the mechanical system in equilibrium. The 
bush can be excited under definite initial conditions, and the energy of the initial excitation 
turns out to be trapped in this bush. 

We found that there exist 18 different by symmetry bushes of vibrational modes for 
the considered mechanical system. We examined the stability of all these bushes for the case 
where interactions between particles are described by the Lennard- Jones potential. The main 
results are listed in Tables [I] and 0. 

It was proved that all of the above mentioned bushes possess stability domains of finite 
size and, therefore, they really can be treated as certain physical objects. 

Appendix 1. The point group Oh and its irreducible 
representations 

The group Oh consists of 48 symmetry elements. We can obtain all these elements with the 
aid of the different products of a certain set of generators. The minimal number of generators 

20 If the centered particle can move, we obtain some additional bushes to those in Tables [ll 0. 



for group Oh is equal 2, but it is more convenient, for our purpose, to use an extended set of 
generators. In notation of the handbook by Kovalev ||, we chose the following five generators: 
h2, h%, /15, hi3, /125. Here h 2 , h% and h\ 3 are 180° rotations about the two-fold axes [100], [010] 
and [110], respectively; /15 is a 240° rotation about the three-fold axis [111], and h 2 ^ is the 
inversion. Using these generators we can obtain, step by step, the following chain of groups: 
£2^2] —* D 2 [h2, h 3 ] — > T[h 2 , h 3 , h 5 ] — > 0[h 2 , h 3 , ^5^13] — > Oh- Each group is obtained from 
the preceding group by adding one new generator.EJ 

The group Oh has ten irreducible representations Tj (j = 1, ...,10). Matrices of the 
generators of this group for all irreps are given below in Table |j. 



Table 4: Matrices of the generators of the group Oh 
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Figure 1: Octahedral mechanical system. 
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